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Recent years have witnessed great interest?^ 
in the quantum spin Hall effect (QSHE) which 
is a new quantum state of matter with nontriv- 
ial topological property due to the scientific im- 
portance as a novel quantum state and the tech- 
nological applications in spintronics. Taking ac- 
count of Si, Ge significant importance as semi- 
conductor material and intense interest in the re- 
alization of QSHE for spintronics, here we inves- 
tigate the spin-orbit opened energy gap and the 
band topology in recently synthesized silicene us- 
ing first-principles calculations. We demonstrate 
that silicene with topologically nontrivial elec- 
tronic structures can realize QSHE by exploiting 
adiabatic continuity and direct calculation of the 
Z2 topological invariant. We predict that QSHE 
in silicene can be observed in an experimentally 
accessible low temperature regime with the spin- 
orbit band gap of 1.55 meV, much higher than 
that of graphene due to large spin-orbit coupling 
and the low-buckled structure. Furthermore, we 
find that the gap will increase to 2.90 meV under 
certain pressure strain. Finally, we also study 
germanium with similar low buckled stable struc- 
ture, and predict that SOC opens a band gap of 
23.9 meV, much higher than the liquid nitrogen 
temperature. 

Quantum spin Hall effect (QSHE) with time reversal 
invariance is gapped in the bulk and conducts charge 
and spin in gapless edge states without dissipation at 
the sample boundaries. The existence of QSHE was first 
proposed in a graphene in which the spin-orbit coupling 
(SOC) opens a band gap at the Dirac point by Kane 
and Mele 1 . But the subsequent work found the SOC 
rather weak, which is in fact a second order process for 
a flat graphene, so the QSHE in a flat graphene can 
only occur at unrealistically low temperature^]. So far, 
there is only one proposal that is able to demonstrate 
QSHE in a real system, which is in two dimensional mer- 
cury telluride-cadmium telluride semiconductor quantum 
wells^ in despite of some theoretic suggestions^. Nev- 
ertheless, both graphene and HgTe quantum wells are 
not good enough to be compatible with the present 
silicon-based electronics industry. As the counterpart of 
graphene^ for silicon, silicene has been shown that a 
low buckled two-dimensional hexagonal structure corre- 
sponds to a stable structure, and there are also evidences 
of graphene-like electronic signature in silicene nanorib- 
bons experimentally^^. Therefore, almost every strik- 
ing exceptional property of graphene could be transferred 
to this innovative material with the extra advantage of 




FIG. 1: The lattice geometry of low buckled silicene. a 

,b, The lattice geometry from side view and top view respec- 
tively Note that A-sublattice (red) and B-sublattice (yellow) 
are not coplanar. c, The first Brillouin zone of silicene and its 
points of high symmetry d, The angle is denned as between 
the Si-Si bond and the Z direction normal to the plane. 

easily being incorporated into the silicon-based micro- 
electronics industry. 

The structure of silicene is shown in Fig. [I] We ob- 
tain the low-buckled geometry of minimum energy and 
stability with lattice constant a = 3.86 A and nearest 
neighbor Si-Si distance d = 2.28 A through structural 
optimization and calculations of phonon spectrum. The 
results agree with the previous workM Compared with 
graphene, the larger Si-Si interatomic distance weakens 
the 7T — 7r overlaps, so it cannot maintain the planar struc- 
ture anymore. This results in a low-buckled structure 
with sp 3 -\ike hybrid orbitals. In Fig. [I] one can define 
the angle between the Si-Si bond and the direction nor- 
mal to the plane. The sp 2 (planar), low-buckled and sp 3 
configurations correspond to 6 = 90, = 101.73 and 
= 109.47 degrees respectively. 

To illustrate the band topology of the low-buckled sil- 
icene, we begin with their graphene analog, planar sil- 
icene, and follow its band structure under an adiabatic 
transformation during which the unstable planar honey- 
comb structure is gradually evolved into the low-buckled 
honeycomb structure. Planar silicene with the same 
structure as graphene should have the similar proper- 
ties. Furthermore, since Si atoms have greater intrinsic 
spin-orbit coupling strength than C atoms, it is natural 
to conceive that the quantum spin Hall effect is more sig- 
nificant in planar silicene. According to symmetry, the 
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low energy effective Hamiltonian with SOC in silicene in 
the vicinity of Dirac point K can be described by 

rr[K] _ ( icr z v F (k x + ik y ) \ f , 

where v F is the Fermi velocity of tt electrons near the 
Dirac points with the almost linear energy dispersion, and 
<j z is Pauli matrix. The effective SOC £ for planar silicene 
has the explicit form £ w 2£g| A e |/(9V^^) with A e being 
the energy difference between the 3s and 3p orbitals and 
£o the intrinsic spin-orbit coupling strength respectively. 
The parameter V spcr corresponds to the a bond formed 
by the 3s and 3p orbits. The effective Hamiltonian near 
Dirac point K* can be obtained by the time-reversal op- 
eration on the one near K. The above equation results 
in a spectrum E(k) = ±y/ {vpk) 2 + £ 2 . Therefore, one 
can estimate the energy gap, which is 2£ at the Dirac 
points, to be about the order of 0.1 meV by taking the 
values of the corresponding parameters 15 . Notice that 
in planar silicene (6 = 90) tt orbitals and a orbitals are 
coupled only through the intrinsic SOC. So, the effective 
SOC is in fact a second order process. However, with the 
deviation (0 > 90) from the planar geometry, tt orbitals 
and a orbitals can also directly hybridize. Consequently, 
the magnitude of the effective SOC depends on the angle 
0. As can be expected with increasing the degree of de- 
viation from the planar structure, the effective SOC will 
be incremental, and QSHE will be more significant. 

The argument above is supported by our first- 
principles calculations based on density-functional the- 
ory (DFT). The relativistic electronic structure of sil- 
icene is obtained self-consistently by using the projec- 
tor augmented wave (PAW) pseudopotential method 
implemented in the VASP package^. The exchange- 
correlation potential is treated by Perdew-Burke- 
Ernzerhof (PBE) potential. 

We carry out detailed and systematic calculations of 
the band structure in adiabatic evolvement from the pla- 
nar honeycomb geometry to the low-buckled honeycomb 
geometry. The evolvement of the gap opened by SOC 
for the tt orbital at the Dirac point K from the planar 
honeycomb geometry to the low-buckled honeycomb ge- 
ometry is shown in Fig. |2^i. Fig. [2Jd and Fig. |2]i show the 
band structures of planar and low-buckled silicene respec- 
tively with the corresponding structures in Fig. [2^i. The 
band structures of planar and low-buckled geometry are 
slightly different in consideration of that the gap induced 
by the effective SOC increases and the degeneracies at 
some k points split. The difference of the both band 
structures along TM high symmetric orientation and in 
the energy range from —3 eV to —2 eV actually means 
that a orbital and tt orbital can directly hybridize only 
in low-buckled geometry. We can find that the gap in- 
duced by SOC for a orbitals is 34.0 meV at T point in 
both geometries. As is also shown in the figure that the 
magnitude of the gap induced by effective SOC for the tt 
orbital at the K point in planar geometry is 0.07 meV, 
which is in agreement with the estimate obtained from 
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FIG. 2: The adiabatic evolvement of the gap, calcu- 
lated relativistic band structure and the topological 
Z2 invariant of silicene. a, The evolvement of the gap 
opened by SOC for the tt orbital at the Dirac point K from 
the planar honeycomb geometry to the low-buckled honey- 
comb geometry with keeping Si-Si bond length constant. The 
top and down abscissas correspond to the difference of Z ab- 
scissa between A-sublattice and B-sublattice and the angle 
aforementioned respectively during evolvement. b and d are 
the relativistic band structures with the corresponding geome- 
tries in a. b, d, Main panel: the relativistic band structure 
of planar silicene and low-buckled silicene respectively. Inset: 
zooming in the energy dispersion near the K point and the 
gap induced by SOC. c, The n-field configuration for silicence. 
The calculated torus in Brillouin zone is spanned by G\ and 
Gi- Note that the two reciprocal lattice vector form a angle 
of 120 degrees. The white and black circles denote n — 1 and 
— 1, respectively, while the blank denotes 0. The Z2 invariant 
is 1 obtained by summing the n-field over half of the torus. 



the tight-binding model discussed above. In low-buckled 
structure, the magnitude of the gap is 1.55 meV, which 
corresponds to 18 K. The top coordinate in Fig. [2^i de- 
notes the distance of two nonequivalent Si atoms within 
a primitive cell in the vertical direction. The figure in- 
dicates that the gap with the magnitude of 0.07 meV 
in the planar structure has been increasing to 1.55 meV 
in the low-buckled structure with energy minimum and 
stability. Most importantly, the gap is not closed. There- 
fore, the low-buckled silicene with energy minimum and 
stability must share the same nontrivial topological prop- 
erties as the planar silicene. Consequently, QSHE can be 
realized in the low-buckled silicene, namely the native 
geometry of silicene. The argument can be confirmed by 
direct calculation of the Z2 topological invariant. 

One can interpret nonzero topological Z2 invariant as 
an obstruction to make the wave functions smoothly de- 
fined over half of the entire Brillouin zone under a cer- 
tain gauge with the time reversal constraintP^ni. The 
band topology can be characterized by the Z2 invari- 
ant. Z2 = 1 characterizes a nontrivial band topology 
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FIG. 3: The gap induced by SOC and the Fermi ve- 
locity of charge carriers vf near the Dirac points are 
calculated from first-principles method under differ- 
ent hydrostatic strain condition. The black circles and 
red diamonds mark the gap and the Fermi velocity vf near 
the Dirac points under different hydrostatic strain A. Inset: 
Energy of unit cell versus different hydrostatic strain condi- 
tion. 



while Z2 = means a trivial band topology. Here we 
follow the method in Ref. [21] to directly perform the 
lattice computation of the Z2 invariants from our first- 
principles methocP^EI] The n— field configuration for 
the low-buckled silicene is shown in Fig. [2j; from first- 
principles calculations. It should be noted that differ- 
ent gauge choices result in different n- field configurations, 
however the sum of the n-field over half of the Brillouin 
zone is gauge invariant module 2, namely Z2 topological 
invariant. As shown in Fig. [2^, low-buckled silicene has 
nontrivial band topology with the topological invariant 
Z2 = 1. Therefore, QSHE can be realized in the low- 
buckled silicene, that is the native geometry of silicene. 

In what follows, we investigate the gap opened by SOC 
at Dirac points related to QSHE and the Fermi velocity 
of charge carriers vp near the Dirac points in a series 
of silicene geometries under hydrostatic strain from first- 
principles method. We find that while the largest pres- 
sure strain can reach —6% without destroying the non- 
trivial topological properties of those systems, the mag- 
nitude of the gap at Dirac points induced by SOC can be 
up to 2.90 meV, which corresponds to 34 K. As shown in 
Fig. [3j the magnitude of the gap at Dirac points induced 
by SOC is incremental with the decrease of hydrostatic 
strain A, which is defined as A = (a — ao)/ao x 100%, 
where ao and a being the lattice constant without and 
with hydrostatic strain respectively. In the pressure 
strain range, the QSHE can be also realized in the system 
and even more pronounced. The figure also indicates the 
greater the angle the greater the gap. In addition, we 
evaluate the Fermi velocity of charge carriers vp near the 
Dirac points under different hydrostatic strain and find 
that the magnitude of the hydrostatic strain does not 



FIG. 4: The calculated relativistic band structure of 
germanium with honeycomb structure, a, b, Main 
panel: the relativistic band structure of Germanium with pla- 
nar and low buckled honeycomb structure respectively. Inset: 
zooming in the energy dispersion near the K point and the 
gap induced by SOC. 



significantly change the carrier Fermi velocity vp. The 
value is slightly less than the typical value in graphene, 
say, 10 6 m/s due to the larger Si-Si atomic distance. 

Recentl y, s everal experiments on silicene have been 
reportecP^^. They have not only proven silicene adopt- 
ing slightly buckled honeycomb geometry and possessing 
the band dispersion with a behavior analogous to the 
Dirac cones of graphene but also synthesized a silicene 
sheet through epitaxial growth. With the advancement 
in experimental techniques, we expect that silicene with 
high-qualify will soon be manufactured. The experimen- 
tal data available can be compared with our theoretical 
prediction then. 

Although Germanium with two dimensional honey- 
comb geometry has not been yet synthesized in exper- 
iments so far, we also conduct a detailed study on ger- 
manium with two dimensional honeycomb structure be- 
cause of its similarity to the other group IVA elements 
in the periodic table, as well as its significant importance 
as semiconductor material. After structural optimization 
and calculations of phonon spectrum, the low-buckled ge- 
ometry of minimum energy and stability with lattice con- 
stant a = 4.02 A and nearest neighbor Ge-Ge distance 
d = 2.42 A is obtained. As shown in Fig. [4^1 and Fig. [4)3, 
Ge with low-buckled honeycomb structure is insulator 
while Ge with planar honeycomb structure is metallic. 
Fig. [4)3 indicates that the magnitude of the gap induced 
by effective SOC for the tt orbital at the K point in low- 
buckled geometry is 23.9 meV corresponding to 277 K 
which is much higher than the liquid nitrogen tempera- 
ture. The direct calculation for topological Z2 invariant 
proves that Ge with low-buckled honeycomb structure 
has nontrivial band topology. Therefore, we predict that 
QSHE will be realized in native germanium with two di- 
mensional low-buckled honeycomb geometry and easily 
observed experimentally once this novel material is syn- 
thesized. 

In conclusion, we have shown both silicene and Ge 
with two dimensional honeycomb geometry have nontriv- 
ial topological properties in their native structure. In ad- 
dition, the QSHE in silicene can be more significant under 
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a range of hydrostatic strain due to the increasing gap 
size. These are confirmed by direct calculations of the 
topological Z<2 invariants from first-principles methods. 
Silicene and Ge with the low buckled honeycomb geom- 
etry have the novel physical properties akin to graphene 
such as the linear energy dispersion at the Fermi level. 
Besides, Silicene and Ge with the low buckled geome- 



try and great SOC can be not only synthesized and pro- 
cessed using mature semiconductor techniques but also 
more easily integrated into the current electronics. All 
of these make silicene and Ge with the low buckled hon- 
eycomb geometry cornucopias of fundamental physics in- 
terests and promising applications. 
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